Asymmetric quintuplet condensation in the frustrated S=l spin dimer compound 
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Ba3Mn2 08 is a spin-dimer compound based on pairs of S = 1, 3d 2 , Mn 5+ ions arranged on a 
triangular lattice. Antiferromagnetic intradimer exchange leads to a singlet ground state in zero 
field, with excited triplet and quintuplet states at higher energy. High field thermodynamic measure- 
ments are used to establish the phase diagram, revealing a substantial asymmetry of the quintuplet 
condensate. This striking effect, all but absent for the triplet condensate, is due to a fundamental 
asymmetry in quantum fluctuations of the paramagnetic phases near the various critical fields. 



Spin dimer compounds provide fertile ground for 
studying novel ordered states in the context of a magnetic 
material [![. Typically comprising S=l/2 ions, antifer- 
romagnetic intradimer exchange leads to a groundstate 
that is a product of singlets. Interaction between dimcrs 
causes the excited triplet states to broaden in energy, 
while application of an external magnetic field leads to 
familiar Zeeman splitting. Above a characteristic critical 
field the minimum in the S z = 1 triplon band crosses 
the singlet in energy, leading to long range order at T=0 
0, [H ■ Depending on the balance of potential and kinetic 
energy of the triplets, the ground state can be variously 
approximated as a Bose-Einstein condensate (canted XY 
antiferromagnet) [1, H ^ fil 0| j a triplet crystal (inhomo- 
geneous Ising order) 8| or even a supersolid [9(. Sig- 
nificantly, quantum fluctuations play a prominent role 
in these systems, providing means to dynamically tune 
quantities such as the quasiparticle effective mass (Io| . 

An intriguing facet of spin dimer compounds, which 
has not been widely studied, is the possibility of higher 
moment ordered states at fields beyond that required to 
saturate the triplet states. Considering the example of 
dimers composed of spin 1 ions, the energy spectrum con- 
sists in zero magnetic field of excited triplet and quintu- 
plet states above the singlet ground state (Fig. [lja)). 
Application of a magnetic field splits both the triplet 
and quintuplet states, and depending on the size of in- 
teractions, can lead at low temperatures to a distinct 
S z = 1 triplet/ S z = 2 quintuplet ordered state at fields 
above those required to the saturate the triplet states. 
The organic biradical magnet F2PNNNO, composed of a 
strongly ferromagnetically linked spin 1/2 pair linked an- 
tiferromagnetically to another pair, acts as a spin 1 dimer 
111 . 1 1 21 ] . Experiments have shown two regions of linearly 
increasing magnetization, first as the S z = 1 triplet states 
and later as the S z = 2 quintuplet states are populated, 
in addition to an ordered phase for the first region. In 
this letter we explore the nature of the quintuplet con- 



densate in the novel spin dimer compound BasMi^Os, 
finding a substantial asymmetry that can be accounted 
for by consideration of quantum fluctuations near the 
critical fields. 

Ba3Mn.2 08 crystallizes in the rhombohedral R3m 
structure, and is comprised of pairs of Mn 5+ 3d 2 5 = 1 
ions arranged vertically on hexagonal layers (Fig. QJb)) 
fl3| . Inelastic neutron scattering (INS) measurements 
revealed the dominant intradimer exchange (Jo = 1.642 
meV), as well as a hierarchy of additional exchange inter- 
actions, including the dominant exchanges linking dimcrs 
within a hexagonal plane (J2 — J3 = 0.1136 meV) and 
the dominant out of plane interactions (Ji = 0.118 
meV and J 4 = 0.037) 0, [H]. Electron paramag- 
netic resonance (EPR) experiments in the diluted com- 
pound BastVi-^Mn^Os (where the V 5+ 3d ion car- 
ries no moment) reveal an almost isotropic g-tensor, with 
g cc = 1.96 and g aa = 1.97, and a single ion uniaxial 
anisotropy characterized by D = 5.81 GHz, correspond- 
ing to 0.024 meV Hi'. Similar measurements for the 
pure compound Ba3Mn 2 0s indicate a zero field splitting 
of the triplet characterized by \D\ = 0.032 meV [17J, the 
difference being due to the presence of additional sym- 
metric anisotropics in the concentrated lattice, in par- 
ticular dipolar coupling between the two ions on each 
dimer. Thermodynamic measurements in fields up to 
30T revealed a complex phase diagram associated with 
the triplet states, consisting of two distinct ordered states 
for fields oriented away from the c axis, attributed to 
the competition between single ion anisotropy and inter- 
dimer exchange in the frustrated lattice [181 ]. Here we 
probe the behavior of single crystals of BasM^Os via 
magnetization, heat capacity and magnetocaloric effect 
measurements in fields large enough to close the gap to 
the 5* 2 = 2 quintuplet states. 

Extraction magnetization measurements were per- 
formed in pulsed magnetic fields up to 60T in a 3 He re- 
frigerator for fields applied perpendicular to the c axis 
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FIG. 1: (Color online) (a) Energy spectrum versus field 
for isolated dimer composed of two spin 1 ions. (b) 
Schematic diagram illustrating the coordination of Mn dimers 
in Ba3Mn208- J0-J4 represent exchange interactions, de- 
scribed in main text, (c) Magnetization curves for BasM^Os 
taken for fields applied perpendicular to the c axis. Successive 
temperature sweeps are offset by 0.2 /is for clarity. 



[19{ . Data were obtained by integrating the field deriva- 
tive of the magnetization and are plotted for increasing 
fields in Figure HJc). The data were cross calibrated with 
low field SQUID measurements to attain absolute values 
of the magnetization. Critical fields, evident as disconti- 
nuities in the slope of the magnetization, were determined 
from peaks in the second derivative of magnetization with 
respect to field. At base temperature, the data show no 
magnetization up to .ff c i=8.73T, followed by linearly in- 
creasing magnetization as the dimers are populated with 
triplet states, and a plateau at 1 \ib between i? c2 =26.46T 
and iJ c ,3=32.42T corresponding to one S z = 1 triplet per 
dimer. There is then a second region of linearly increas- 
ing magnetization as the quintuplet band is filled, before 
finally reaching the full saturation magnetization of 2 \ib 
at iJ c4 =47.9T [H. 

The derivative of magnetization with respect to field 
is plotted in Figure [2] for various temperatures. The 
data were obtained by measuring the time derivative of 
both the magnetization and field, and taking the quo- 
tient. The signal to noise ratio diminishes at high fields, 
for which the field derivative tends to 0. These curves, 
shown for increasing fields, reveal clear peaks associated 
with each of the critical fields. A representative curve for 
decreasing field at base temperature is plotted revealing 
the same qualitative behavior. For falling fields the field 
derivative is smaller, resulting in a worse signal to noise 
ratio. The data also show a second peak to the right 



(left) of H c i and H c % (H C 2 and H C 4). These peaks do not 
mark phase transitions but rather signify a zero temper- 
ature crossover between two distinct phases observed in 
our earlier work [l8| . 

Inspection of Fig. [2] reveals that at lowest temperature 
the peak at H c i is significantly higher than the peaks 
at the other three critical fields. The magnitude of the 
peaks at the critical fields is set by the number of low- 
energy excited states available at each critical point. At 
each critical field H cu (where v = {1,2,3,4}), the sys- 
tem can be described as a dilute gas of gaplcss bosons 
with a dispersion relation w v (q) that is quadratic at low 
energies. Since each boson is a quasiparticle that car- 
ries a unit of magnetization (S z = 1) along the field 
direction, the boson density is the magnetization per site 
(Sj). The diagonal components of the mass tensor are 
m~i = 9 2 w !/ (q)/5 2 3 i | q= Q, where i = {x,y, z} and Q is 
the wavevector where lo v (q) is minimized. The enhanced 
peak at H c ± implies a larger density of states, and by 
extension a larger effective mass for the corresponding 
bosons. The difference in the effective mass is a direct 
result of the strong asymmetry in the zero point fluctua- 
tions of the paramagnetic phases. Conservation of total 
S z implies that such fluctuations are absent above H c ±. 
In contrast, zero-point fluctuations consisting of creation 
(annihilation) of S z — 2 quintuplct-singlct pairs (triplet- 
triplet pairs) are present for H C 2 < H < H c ^. Their 
effect is to reduce the mass of the bosons. This can be 
easily appreciated by a thought experiment in which the 
field is fixed at some arbitrary value between H C 2 and 
H c z and the interdimer exchange couplings are increased 
by applying pressure. The applied pressure P will close 
the gap at a critical point (P = P c ) where H C 2 becomes 
equal to H c ^, causing the single-particle excitation spec- 
trum to become linear [z = 1). Such an XY quantum 
phase transition in O = rf + 1 = 4 dimensions is driven 
by the zero-point (phase) fluctuations under considera- 
tion 2l|. A simple continuity argument implies that the 



boson mass has to be reduced by such fluctuations in 
order to get a massless spectrum when they diverge at 

P = Pc- 

In the dilute limit relevant for the regions around each 
of the quantum critical points, the inter-particle dis- 
tance p -1 / 3 (where p is the boson density) is much larger 
than the scattering length a (lattice parameter) of the 
hard-core repulsion between bosons. The effective boson- 
boson repulsion, vp = To(Q, Q), is obtained from a lad- 
der summation 12211: 
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p)+u v (k' -p) 
(1) 

(J-2 + J 3 ) 7q + ir7 q + if 7 q , where U -> 00 
comes from the hard-core repulsion, the rest of the terms 
come from nearest neighbor repulsion terms and 7 q 's are 



Here Vq = U 
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FIG. 2: (Color online) Field derivative of magnetization of 
Ba3Mn2 08 for increasing fields applied perpendicular to the 
c axis. Successive temperature sweeps offset by 0.0045 /Lts/T 
for clarity. Dotted line shows representative magnetization 
derivative for decreasing fields at 0.345K, revealing no signif- 
icant hysteresis. 

defined as follows: 

7q = cos(g 3 ) + cos(q 3 - q{) + cos (g 3 - q 2 ) 

7q = cos(gi) + cos(r/ 2 ) + cos (gi - q 2 ) 

7q = cos (q 3 +q\- q 2 ) + cos (q 3 -qi + q 2 ) 

+cos (q 3 - qx - q 2 ) , (2) 

where q n 's refer to the principal reciprocal axes. The 
dispersions w"(q) were computed by using a generalized 
spin- wave approach [23[ . For the paramagnetic (PM) re- 
gion near each QCP, a mean field treatment of the effec- 
tive interaction vq leads to a rcnormalizcd chemical po- 
tential fi v = jtx„ — 2t>op with /Ay = gfiB (— 1)" {H cv — H) 
Low temperature thermodynamic properties can be 
computed from the renormalized chemical potential in 
the PM region around each QCP. The striking asym- 
metry between H C 4 and the rest of the critical fields 
follows from the fact that the corresponding masses 
m 4 || = 1.9357/v" 1 and ni4± = 0.4531X -1 arc roughly 
two times bigger than the masses at the other three crit- 
ical fields: m 3 || = 0.7834X" 1 , m 3± = 0.1833X m 2 \\ = 
QmilK' 1 , m 2A _ = 0.1633X" 1 , and m n = 0.9133R- 1 , 
mu = 0.2138R- 1 . 

In general, the mass asymmetry present in BasM^Os 
induces asymmetries in all of the thermodynamic proper- 
ties. For instance, the condensation temperature (/i(T = 
T c ) = 0) near the critical fields is: T c oc p 2 l Z Trv u 1 oc 
\H — Hcv^^mZ 1 ^ 3 , where we have used that Tq(Q, Q) oc 
m~ 1 . The increased effective mass at H C 4 reduces the or- 
dering temperature, causing the phase boundary to be 
steeper near H c3 than near H c ±. In general, all spin 
dimer compounds should exhibit some degree of asym- 
metry because zero-point fluctuations arc absent only in 
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FIG. 3: (Color online) (a) MCE curves plotted versus field 
for increasing field (solid black lines) and decreasing fields 
(dashed red lines) for fields applied perpendicular to the c 
axis. Phase transitions can be determined from peaks in 
derivative with respect to field, and are shown as open trian- 
gles. Vertical shading indicates broad maximum associated 
with rapidly changing temperature derivative of magnetiza- 
tion at the critical field, (b) Heat capacity, taken on cooling, 
as a function of temperature for fields applied perpendicular 
to the c axis. 



the highest field paramagnetic phase. The effect is espe- 
cially striking in Ba 3 Mn208 because H c3 — H c2 ~ 6T is 
considerably smaller than H c i — H c3 ~ 15. 5T, implying 
the gap in energy to create a zero-point fluctuation near 
H c3 is relatively small. 

To observe the additional consequences of the mass 
asymmetry, we measured the high field phase diagram 
by both magnetocaloric effect (MCE) and heat capacity 
experiments performed in a 45T hybrid superconduct- 
ing/resistive magnet at 3 Hc temperatures for fields per- 
pendicular to the c axis. For MCE measurements, tran- 
sitions are marked by a sharp increase (decrease) in tem- 
perature on entering (leaving) the ordered phase, while 
in practice the phase boundary is determined from peaks 
in the derivative (see inset to Fig. [3)Ja)). The verti- 
cal shaded region in Fig. [Ha) refers to a temperature 
independent effect, associated with a rapidly changing 
temperature derivative of magnetization at H c3 (l8j . 

Heat capacity measurements were performed using 
standard thermal relaxation time calorimctry, in the 
large AT limit. Accordingly, C p is calculated from the 
instantaneous time derivative of the calorimeter tempera- 
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FIG. 4: (Color online) Phase diagram for fields applied per- 
pendicular to the c axis. Representative magnetization curve 
for T = 0.50 K is shown for reference (right axis). Blue re- 
gion corresponds to the quintuplet condensate. Solid and 
open circles correspond to heat capacity data taken on heating 
and cooling respectively, illustrating the apparent hysteresis 
caused by the experimental technique. Dark and light yel- 
low regions correspond to phase II and phase I of the triplet 
condensate as defined in ref. Hal. 



ture and the conductance of the thermal link as discussed 
elsewhere [2~4| . For these experiments, the primary ben- 
efit of this technique is its expediency, which unfortu- 
nately comes at the cost of reduced accuracy. The re- 
duced accuracy caused a small extraneous hysteresis on 
heating versus cooling. To improve the accuracy, multi- 
ple temperature curves were measured at each field and 
then averaged. The data taken on cooling are plotted in 
Fig. [3] They reveal an increasing critical temperature 
for fields leading up to 37T which has the maximum T c 
of 0.63K. Above 37T the critical temperature begins to 
decrease again. T c values extracted from heat capacity 
data for both cooling and heating are shown as a lower 
and upper bound on the phase boundary in Fig. [4] The 
integrated entropy under the curve appears to track the 
critical temperature, so that a higher T c implies more 
entropy This effect is also produced by the asymmetry 
in the zero-point fluctuations that shift the maximum of 
the order parameter towards the critical field with larg er 
fluctuations (i? C 3 for the triplet-quintuplet phase) [23| . 

The entire phase diagram for fields applied perpendic- 
ular to the c axis is shown in Figure [4] together with a 
magnetization curve, taken at 0.50K, showing the four 
critical fields as discontinuities in the slope. All told 
the phase diagram reveals at least three distinct ordered 
states across a large field range. Inspection of this figure 
reveals a striking asymmetry in the quintuplet conden- 
sate, for which the maximum T c of the quinton phase is 
found at -ff=37T, less than 5T from H C 3 but more than 
10T from H C 4. Additionally, a small degree of asym- 
metry was found in the ordered states of the singlet- 
triplet regime, with phase II slightly narrower on the 
high field side than the low field side. Both asymmetries 
are expected from the different masses reported above: 
vniijmu ~ 0.7639 and ma/msi ~ 2.4711. The measured 
ratio between maximum T c 's, T£* ax /T^ ax ~ 1.37, agrees 
well with the ratio of 4/3 between the effective intcrdimcr 



exchange interactions obtained from the two level model 
18 1 for the singlet-triplet and triplet-quintuplet phases 
23|. 

In summary, we have determined the high field phase 
diagram of the S = 1 dimer compound BaaMi^Og. 
Quantum fluctuations play a profound role in de- 
termining the relative symmetry/asymmetry of the 
triplet/quintuplet condensates, respectively. This mate- 
rial provides an ideal test bed for studying and tuning 
the effect of zero-point fluctuations. 
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